Critical Care

OBJECTIVE: Sepsis-associated encephalopathy (SAE) and sepsis-associated delirium (SAD) lack confirmatory laboratory tests or biomarkers. We evaluated whether the Th17/interleukin (IL)-17 axis or neuroglial-injury biomarkers reflect acute brain dysfunction beyond systemic inflammation and whether biomarker profiles differ between these phenotypes. DATA SOURCES: EMBASE, MEDLINE/PubMed, and the Cochrane Library from inception to August 1, 2026, supplemented by backward and forward citation searching, following PRISMA 2020. STUDY SELECTION: Observational studies of adults with sepsis classified by SAE or SAD status and reporting circulating biomarker concentrations by neurological group or biomarker-outcome associations. Eighteen studies (2,388 participants; 1,863 with sepsis) were included. DATA EXTRACTION: Two reviewers independently extracted study, phenotype, assay, and outcome data. Risk of bias was assessed with the Newcastle-Ottawa Scale and certainty with GRADE. DATA SYNTHESIS: Random-effects models pooled standardized mean differences (SMDs) separately against non-septic and septic controls and odds ratios (ORs) for biomarker-outcome associations. IL-6 was the only mediator elevated in both comparative frameworks and associated with outcomes (non-septic controls: SMD, 0.54 [0.25-0.82]; septic controls: SMD, 1.05 [0.51-1.59]; OR, 1.85 [1.50-2.29]). S100B, neuron-specific enolase, IL-10, tumor necrosis factor; and interferon-gamma; were elevated versus septic controls, whereas IL-17 was modestly elevated only versus non-septic controls (SMD, 0.37 [0.01-0.73]). The IL-6 septic-control effect was concentrated in SAE (p<0.0001) and remained robust in leave-one-out and risk-of-bias sensitivity analyses. Certainty was low to very low, and all 14 computable prediction intervals included the null. CONCLUSIONS: IL-6 emerges as the leading circulating candidate biomarker of sepsis-associated brain dysfunction, while S100B and NSE capture neuroglial injury beyond sepsis. Current evidence does not establish a discriminating role for the Th17 axis. Future studies require standardized neurological phenotyping and appropriately matched septic controls to establish diagnostic value.

BACKGROUND Traditional organ dysfunction classifications after severe trauma rely on binary categories (e.g., chronic critical illness [CCI] versus rapid recovery), which likely oversimplify the true heterogeneity of postinjury organ dysfunction trajectories and their underlying immune-biologic drivers. We hypothesized that severely injured trauma patients demonstrate multiple, distinct organ dysfunction phenotypes distinguished by unique immune signatures. METHODS We performed a prospective, longitudinal, observational cohort study of severely injured trauma patients at high risk for multiple organ dysfunction at a single Level 1 trauma center. Organ dysfunction was quantified daily through 28 days, death, or discharge using a trauma-specific modified Sequential Organ Failure Assessment (mSOFA) score excluding the neurologic subdomain to isolate systemic organ dysfunction from traumatic brain injury. We utilized group-based trajectory modeling to identify longitudinal mSOFA phenotypes. Plasma immune protein trajectories were compared across phenotypes using mixed-effects models with time interaction. RESULTS Sixty-nine patients were enrolled. Trajectory modeling identified five mSOFA phenotypes (c1-c5) with differences in resuscitation requirements, shock severity, and 28-day mortality, but not baseline demographics, injury severity, or early physiologic parameters. The phenotypes revealed organ dysfunction heterogeneity not fully captured by CCI/rapid recovery classification. Forty of 92 measured immune markers showed phenotype-differentiated trajectories (FDR<0.1). The most fulminant phenotype (c5) was distinguished by IL-17C, CXCL5, and CXCL6, implicating epithelial-barrier and neutrophil-trafficking mechanisms. The phenotype with highest CCI incidence (c3) was distinguished by rising CXCL9, ARTN, and TWEAK, implicating Th1 polarization and tissue-inflammation/repair signaling. The clinically distinct, older-age phenotype (c4) carried the largest differentiating signature (29 proteins), spanning renal injury/repair, endothelial activation, and adaptive lymphocyte modulation. CONCLUSIONS Organ dysfunction after severe trauma is more biologically heterogeneous than traditional clinical classifications suggest. Distinct mSOFA trajectory phenotypes correspond to distinct immune signatures, supporting separable pathobiology rather than a single axis of injury severity, with implications for phenotype-targeted precision immunomodulatory approaches in critical illness.

Incidence of ventilator-associated pneumonia (VAP) is increased in patients admitted for brain injury. Staphylococcus aureus appears to be a particularly common pathogen in these respiratory infections. Our aim was to confirm the extent of S. aureus involvement in VAP among patients with neurological injuries. We conducted a retrospective, multicenter study including patients with and without neurological injuries who developed VAP during their ICU stay. Among 15 702 patients admitted in the participant centers during the study period. 407 were finally included in the study. S. aureus was involved in 122 (61.6%) VAP in brain-injured patients and 63 (30.1%) of non-brain-injured patients. In brain-injured patients, multivariate analysis identified age (0.97 [0.95; 0.98]), antibioprophylaxis (0.31 [0.18; 0.53]; p < 0.0001), subarachnoid hemorrhage (2.18 [1.01; 4.68]; p=0.0461), lower Glasgow coma score (0.92 [0.87; 0.97]; p= 0.004) and SAPS II (0.98 [0.97; 1]; p=0.0174) as critical factors associated with the risk of S. aureus in VAP. In the non-brain-injured population, multivariate analysis identified antibioprophylaxis (0.38 [0.22; 0.65]; p=0.0006), age (0.96 [0.95; 0.98]; p<0.0001), SAPS II (0.98 [0.97; 1]; p=0.0265), thoracic surgery (0.21 [0.09; 0.48]; p=0.0002) and peripheral neurosurgery (1.95 [1.2; 3.19]; p=0.0076) as critical factors. A total of 172 patients (42.3%) died during their ICU stay and 90-days mortality was 194 (47.7%). No difference of mortality was noticed function of VAP-associated pathogen.

Introduction: Fighting the ventilator (FV) is characterized by a patients active counteraction of respiratory gas flow across multiple breaths during mechanical ventilation. An effective ventilation is prevented during such episodes. The objective of this exploratory study was to evaluate clinicians knowledge, perception and agreement regarding this condition. Methods: In a cross-sectional online survey, medical professionals with experience in intensive care medicine (mostly neonatal and pediatric patients) responded to 11 questions covering various clinical aspects of FV. Results: All respondents were familiar with the phenomenon and assessed it to be clinically relevant. Neonates were perceived to be predominantly affected. Promoting factors of FV included inadequate sedation and ventilation-associated factors. In clinical practice, the participants identified FV by clinical evaluation, analysis of ventilator waveforms and capnography and used both pharmacological and non-pharmacological intervention strategies. There was no consensus concerning the English terminology of the condition. Only few respondents were familiar with local guidelines to manage FV. Conclusion: Although there was consensus on the relevance of FV, diagnostics and management strategies were highly heterogeneous. Many participants agreed that poor patient-ventilator synchrony was an important cause of FV. Future work should focus on characterizing and generating evidence on FV.

Introduction and objectives Hydrocephalus is a common complication of aneurysmal subarachnoid haemorrhage. This is one of several standard indications for the insertion of an external ventricular drain, but there are no standard guidelines for timing of extraventricular drain discontinuation in this population. Therefore, the study aims were to examine the 1) factors associated with extraventricular drain duration in aneurysmal subarachnoid haemorrhage and 2) association between extraventricular drain duration and clinical and resource-based outcomes. Methods In this single centre retrospective cohort study, data for adult patients treated in the ICU between 2015-2023 with aneurysmal subarachnoid haemorrhage and requiring an extraventricular drain were analysed. Data were extracted from both a prospective patient registry and the electronic medical record. Linear regression modelling was used to identify the factors associated with extraventricular drain duration. Measures used to examine the association between extraventricular drain duration, and the outcomes of interest included hazard ratio, relative risk, and mean differences as appropriate for a specific outcome. Results Of the participants (n=205), 68% (n = 140) were female, mean age was 60.7 (SD 13.3) years and 71% (n = 145) had a modified Fisher score of 4. Mean extraventricular drain duration was 15.3 (SD 8.6) days. The number of days mechanically ventilated and extraventricular drain challenge by clamp trial were associated with an increased duration of extraventricular drain treatment. Increased extraventricular drain duration was associated with an increased incidence of ventriculoperitoneal shunt insertion, and an increase in ICU and hospital length of stay when extraventricular drain duration exceeded 20 days. Conclusion The factors associated with an increased duration of extraventricular drain treatment were a longer duration of mechanical ventilation. The outcomes associated with an increased duration of extraventricular drain treatment were an increased incidence of ventriculoperitoneal shunt insertion, an increased incidence of ventriculitis when extraventricular drain duration exceeded 20 days, and an increase in ICU length of stay. These study findings may inform the conduct of a larger multi-centre study to examine the factors associated with duration of treatment with extraventricular drain, and outcomes associated with extraventricular drain duration.

Objectives: Using a nationally representative database of U.S. hospitalizations, we aimed to (1) characterize contemporary trends in hypoxic-ischemic encephalopathy (HIE) prevalence by severity, therapeutic hypothermia (TH) utilization, and associated in-hospital mortality; and (2) examine the association between TH and mortality by HIE severity and gestational age (GA), focusing on 35- and 36-weeks GA. Methods: We analyzed newborn hospitalizations from the National Inpatient Sample, 2016-2023. Unspecified, mild, moderate, and severe HIE were identified using ICD-10-CM diagnosis codes P91.60-P91.63; and 35- and 36-weeks GA using codes P07.38-P07.39. TH was identified using ICD-10-PCS codes 6A4Z0ZZ and 6A4Z1ZZ. Survey-weighted multivariable generalized linear models with marginal standardization were used to estimate adjusted mortality risk differences (RD) with surprisal (S) values, and assess whether associations between TH and mortality differed by HIE severity across GA strata. Results: Among 29,767,238 hospital discharges, HIE was identified in 66,535 neonates (2.24 per 1,000), and TH used in 19.2% of HIE cases. In-hospital mortality was 9.74% among all HIE cases and 6.65% among TH-treated neonates. From 2016-23, HIE prevalence nearly doubled (1.58 to 3.01 per 1,000; p-for-trend <0.001) and TH utilization more than doubled (0.33 to 0.74 per 1,000; p-for-trend <0.001), while TH use as a proportion of HIE remained stable (17.8%-20.8%; p-for-trend=0.39). HIE severity distribution was: mild (15%), moderate (29.4%), severe (18.3%), and unspecified (37.3%), with moderate HIE increasing over time (22.3%-33.5%). Crude mortality was consistently highest for severe HIE (range: 33.65%-35.25%). Among infants with severe HIE, TH was associated with reduced mortality overall (RD -10.12%; 95% CI -14.86, -5.39; p<0.001; S=15.12), and specifically at 36-weeks GA (RD -21.98%; 95% CI -35.85, -8.12; p=0.002; S=9.06), but not at 35-weeks (RD -1.83%; 95% CI -21.2, 17.55; p=0.86; S=0.22). Conclusions: TH was associated with mortality reduction in severe HIE overall and at 36-weeks GA, but this benefit remains uncertain at 35-weeks GA.

Sepsis is a global health problem. Patients in the Intensive Care Unit (ICU) are at risk of having sepsis because of their vulnerability and use of invasive medical devices. Emerging antimicrobial-resistance complicates treatment. Understanding bacterial profiles and antibiotic sensitivity patterns in the ICU plays an important role in managing sepsis, optimizing antimicrobial use, and improving patient outcomes. The aim of this study was to show bacterial profiles and antibiotic sensitivity patterns among patients who were admitted to the ICU at Mbeya Zonal Referral Hospital for 12 months. This was a prospective cross-sectional study among 187 sepsis-suspected patients. Data analysis was performed using SPSS Version 27. Categorical data were summarized by frequency and percentage. Fishers exact test, chi-square test, and modified Poisson regression were used to show the association between patients clinical variables and culture positivity at a 95% CI. Statistical significance was set at a p-value < 0.05. Out of 187 samples, 41.7% were culture positive. Gram-negative bacteria were commonly isolated (65.4%, with Klebsiella pneumoniae and Escherichia coli being commonly isolated. Commonly isolated Gram-positive bacteria were Coagulase-negative Staphylococcus species (19%). 100% of gram-positive bacteria and 96% of gram-negative bacteria were Multi-Drug Resistant (MDR), respectively. In multivariable analysis, patients who were intubated and mechanically ventilated (aPR = 2.208, 95% CI: 1.133-4.305, p = 0.020) and post-laparotomy patients (aPR = 3.241,95% CI: 1.320-7.958, p = 0.010) independently had a higher proportion of Gram-negative bacteria. The study reveals a high level of resistance to commonly used antimicrobials in our setting, necessitating strengthening of the antimicrobial stewardship program and routine surveillance of antimicrobial resistance patterns.

Polyendocrine Metabolic Ovarian Syndrome (PMOS, formerly PCOS) is associated with chronic low-grade inflammation. This study aimed to characterize the peripheral cytokine landscape in women with PMOS and identify cytokines associated with disease pathogenesis and phenotypic variation. Plasma concentrations of 48 cytokines were quantified using a multiplex bead-based assay in 102 women with PMOS and 114 healthy controls. We performed receiver operating characteristic (ROC) analysis, odds ratio estimation, correlation analysis, and network analyses to characterise cytokine alterations and identify phenotype-associated inflammatory signatures. Fourteen cytokines were significantly elevated in PMOS compared with healthy controls: IL-18, IL-1Ra, IL-12p70, IL-13, IL-16, MCP-1, MCP-3, CTACK, Eotaxin, RANTES, IP-10, M-CSF, HGF, and {beta}-NGF. ROC analysis identified IL-16, CTACK, and IP-10 as potential discriminators of phenotype A, whereas IL-18, M-CSF, and HGF demonstrated promising discriminatory performance for phenotype B. Odds ratio analysis identified 11 cytokines with odds ratios greater than 1, with IL-16 showing the strongest association with PMOS. Correlation analysis revealed 10 strongly correlated cytokine pairs (r [&ge;] 0.7) in the PMOS group, suggesting coordinated regulation of inflammatory mediators. Gene ontology and network analysis highlighted alterations in chemotaxis and immune cell migration pathways, implicating these processes in PMOS-associated immune dysregulation. These findings reveal a distinct peripheral inflammatory cytokine signature and phenotype-associated differences in PMOS. IL-16 emerged as a prominent candidate for further investigation, given its strong association with PMOS. Peripheral cytokine profiling may provide a minimally invasive approach to characterising immune alterations and identifying candidate biomarkers for disease characterisation and phenotype stratification.

Abstract Background: Placental insufficiency (PI), characterized by impaired placental function, is a leading cause of stillbirth. Although rare coding variants are known to cause a subset of stillbirths, their contribution to PI is poorly understood. Methods: We conducted a secondary analysis of stillbirth cases (N=292) of the Stillbirth Collaborative Research Network study with exome sequencing available on postmortem placental or fetal samples. PI was defined by standardized placental pathological lesions (e.g., maternal vascular malperfusion) and prespecified clinical criteria for placental dysfunction. Candidate gene-based rare-variant burden analysis was performed using the Sequence Kernel Association Test-Optimal framework, focusing on evolutionarily constrained developmental and placental-function genes. Primary analysis compared stillbirths with and without placental vascular malperfusion. Similar comparisons by clinical PI and other placental pathological lesions were conducted in secondary analysis. Results: Stillbirths with PI phenotypes demonstrated different patterns of genetic variation than stillbirths without PI. Specifically, we observed nominal differential enrichment of constrained variants in biologically prioritized placental-function genes, regulating extravillous trophoblast invasion and spiral artery remodeling (CTNND2, PHACTR1, SLIT1, VEGFA), placental angiogenesis (IGF2R, TEK), mitochondrial homeostasis (OPA1), and fetal growth (DLK1, GNAS). Secondary analyses on specific placental pathological lesions suggested recurrent enrichment of genes within similar placental developmental pathways. Conclusions: We observed nominal differential enrichment of constrained developmental and placental-function genes in stillbirths with PI compared to stillbirths without PI. The observed genetic signals across PI phenotypes suggest shared biological pathways underlie placental dysfunction. Our findings advance genotype-phenotype characterization and help prioritize candidate pathways for larger collaborative genomic and functional studies to improve stillbirth risk stratification and prevention.

Mitochondrial disease arises from variants in either the mitochondrial or nuclear genome that impair mitochondrial oxidative phosphorylation. Complex III (ubiquinol:cytochrome c oxidoreductase, CIII2) is a key enzyme of the mitochondrial respiratory chain, catalyzing electron transfer from reduced coenzyme Q10 (ubiquinol) to complex IV, via cytochrome c. CIII2 has one subunit encoded by mitochondrial DNA, the remainder being nuclear encoded. Individuals with CIII2 deficiency present with significant clinical heterogeneity, on account of pathogenic variants in genes encoding structural CIII2 subunits or assembly factors. UQCC1 (GenBank: NM_018244.5) encodes the mitochondrial complex III assembly factor Ubiquinol-cytochrome c reductase complex assembly factor 1, which is currently not associated with human disease. Here we describe nine subjects from eight unrelated families harboring biallelic UQCC1 variants with heterogeneous, early-onset mitochondrial disease presentations and lactic acidosis. Subject-derived fibroblasts demonstrated a decrease in the steady-state levels of UQCC1 and its interacting partner UQCC2, with structural modelling supportive of variant pathogenicity. BN-PAGE and quantitative proteomics demonstrated decreased CIII2 abundance whilst mass spectrometry (MS)-complexome profiling corroborated these findings, identifying impaired respiratory supercomplex formation and a delay in the incorporation of newly synthesized CIII2 subunits into native assemblies. Lentiviral transduction experiments in subject-derived fibroblasts restored the levels of both UQCC1 protein and CIII2 activity. Collectively, these findings provide convergent genetic, biochemical, proteomic and functional evidence that bi-allelic UQCC1 variants disrupt CIII2 assembly, establishing UQCC1 deficiency as a cause of mitochondrial CIII2 deficiency.

Giant cell arteritis (GCA) is an immune-mediated vasculitis that can cause serious morbidity including sight loss, stroke and aortic aneurysm. There is an unmet need for blood-based biomarkers and novel therapeutics. We measured plasma protein levels using the SomaScan v4.1 platform, comprising 7,273 protein measurements (SOMAmers) reflecting 6,596 unique protein targets. We profiled the proteome in 113 patients who presented acutely with possible GCA. 56 had biopsy-proven GCA and 57 did not have GCA (non-GCA). We found 144 SOMAmers increased and 151 reduced in GCA. Dysregulated biological pathways included IL6-JAK-STAT3 signalling, chemokines for both myeloid (e.g. CXCL2, CXCL3, CCL18, CCL23) and lymphoid cells (e.g. CXCL9, CXCL13), and the complement system (e.g. C9, CFHR5, C5a, Factor B). Enzyme-Linked Immunosorbent Assays (ELISAs) confirmed elevated plasma levels of C5a, the terminal complement complex, Factor B and CFHR5 in GCA. We developed a 26-protein signature that distinguished GCA from non-GCA (area-under-the curve, AUC, 0.93) and validated the diagnostic signature in external datasets. Our data implicate terminal complement pathway activation in GCA, raising the possibility of complement-directed therapies. We show proof-of-principle that a multiplex protein biomarker could be used as a diagnostic adjunct in suspected GCA.

Objective. Temperature readings are routinely stratified by race and ethnicity, though significant group differences may reflect physiology, instrument noise, or categorization rather than biology. This study tests whether these categories partition facial infrared thermography measurements beyond individual variation and baseline device noise. Approach. A secondary analysis of thermography data (n = 951 adults) compares inner-canthus maximum (TCEmax) and full-face maximum (Tmax) temperatures across five administrative categories via one-way analysis of variance (ANOVA), mixed-effects modeling, Bland-Altman analysis against oral reference, and Monte Carlo sensitivity analysis. Main results. ANOVA returns significant group differences (p < 0.001; eta-squared = 0.04 for TCEmax, 0.03 for Tmax), but the total spread (0.2 to 0.3 degrees Celsius) is comparable to device accuracy and within-subject variability. Individual variation accounts for 95% of TCEmax variance, and adding simulated Gaussian noise to participant-level temperatures (sigma = 0.10 to 0.30 degrees Celsius) preserves the complete group-mean rank order in as few as 37.5% of iterations. Significance. Testing whether a group difference survives perturbation at the scale of instrument accuracy offers a check against treating significance alone as evidence of physiology. Here, observed group differences are small relative to measurement variability, and their physiological interpretation remains uncertain, with implications for clinical tools that treat administrative categories as biologically meaningful.

Objective: To determine whether plasma neurofilament light chain (NfL) and glial fibrillary acidic protein (GFAP) are elevated in hospitalised patients with delirium relative to patients with dementia and primary psychiatric disorders (PPD). Methods and analysis: This was a cross-sectional observational study conducted at the Royal Melbourne Hospital, Melbourne, Australia, between April and October 2024, with comparator samples from a research biobank. In Phase 1 (Delirium Discovery Cohort), we measured plasma NfL and GFAP using the Single molecule array (Simoa) platform in hospitalised patients with delirium and compared them with dementia and PPD groups. Biomarker concentrations were the primary outcome and were compared using general linear models adjusted for age, sex, and estimated glomerular filtration rate, with bootstrapped confidence intervals. Discriminative performance was assessed as a secondary outcome using receiver operating characteristic (ROC) curve analysis. In Phase 2, we applied this threshold to independent dementia and PPD groups. Results: In Phase 1 (n =176), inpatients with delirium showed higher concentrations of both NfL and GFAP than patients with dementia ({beta} = 0.78; 95% CI [0.37, 1.20] and {beta} = 0.39 [0.07, 0.73] respectively) or PPD groups ({beta} = 1.17; [0.69, 1.66] and {beta} = 0.83 [0.51, 1.24] respectively). ROC curve analysis for differentiating delirium from dementia and PPD groups yielded an optimal NfL threshold of 24.1 pg/mL (AUC = 0.89, 96% sensitivity, 70% specificity). In comparison, GFAP showed lower accuracy (AUC = 0.85), and combining NfL with GFAP (AUC = 0.90) did not significantly improve performance. In Phase 2 (n = 154), the NfL threshold of > 24 pg/ml maintained strong performance (AUC = 0.95, 96% sensitivity, and 77% specificity). Conclusions: Plasma NfL shows promise as a biomarker for delirium in hospitalised patients, distinguishing it with high sensitivity from patients with dementia and PPD. Larger multicentre studies are needed to refine cut-offs, evaluate temporal dynamics, and determine integration into acute hospital delirium pathways.

Objective A diversity of approaches and cohorts have been published in research literature for suicide risk estimation, leading to some uncertainty around which factors and conditions are most strongly contributing to risk. Methods exist for weighting estimates by propensity score, removing simplifying the interpretation of model effects. Few studies have explored the use of these tools to assess when the effect of certain factors is uncertain or inconsistent. Method A mix of 1,001,737 suicide decedents and living patients were evaluated for suicide risk, relying on a subset of commonly identified factors: demographics, psychiatric conditions, substance use disorders, and social needs. Care utilization patterns were then assessed controlling for these factors in a hospital care setting. Four methods were considered for controlling multivariate confounding of odds predicting death by suicide: multivariate logistic regression, 1-to-1 sample matching, inverse propensity weighting (IPTW), and entropy balancing (EBAL). Results Test statistics showed suicide decedents were much more likely to have various psychiatric and substance use diagnoses, social needs and more likely to heightened care utilization across sampling intervals. Greater utilization was found close to date of death. Effect sizes after IPTW and EBAL weighting were also generally greater than would have been documented by regression alone. Conclusions Propensity weighting is a useful alternative to exploring nuanced effects, especially when the direction and strength of an association differ from traditional statistics. By weighting, the overall imbalance of cases remains consistent and effect terms might reflect more nuanced associations than regression or sample matching alone.

Host genetic variation influences heterogeneity in infectious disease severity. Genome-wide association studies (GWAS) have identified 69 lead single nucleotide polymorphisms (SNPs) associated with critical COVID-19, but most lack known regulatory mechanisms. To date, the genetic cis-regulatory effects of critical COVID-19 GWAS lead SNPs have not been systematically characterized across acute illness and convalescent phases in a COVID-19 specific cohort. Here, performing severity-stratified, longitudinal cis-expression quantitative trait loci (cis-eQTL) mapping and colocalization with critical COVID-19 GWAS signals, we find that host cis-regulatory effects vary across clinical severity and disease phase. We analyzed peripheral blood mononuclear cells from 1,031 patients in the Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC), a prospective multicenter study of adults hospitalized with COVID-19, profiled across 10 timepoints spanning acute illness and convalescence. We identified 23 GWAS lead SNPs that colocalized with cis-eQTLs at 26 candidate causal genes, with 16 detectable only when severity stratification and longitudinal sampling were incorporated. A colocalized, protective allele (chr19:4,717,660) at the TMIGD2 locus upregulated TMIGD2 expression exclusively during the acute phase in the prolonged-critical/fatal cohort, with no detectable effect in the rapid-recovery cohort. Our results show that host genetic effects on critical COVID-19 can be context-dependent, with cis-regulatory effects varying across clinical severity and disease phase. Severity-stratified, longitudinal cis-eQTL designs have greater sensitivity to identify candidate causal genes for critical COVID-19 GWAS associations than single-timepoint cross-sectional studies.

Glycemic control is key to prevention of complications of diabetes. Although eosinophils have been implicated in metabolic regulation, most human studies have been cross-sectional and have yielded inconsistent findings. We conducted a retrospective cohort study of adults with uncontrolled diabetes receiving primary care at Mass General Brigham. We evaluated the relationship between eosinophil count and HbA1c over one year, adjusting for demographic, socioeconomic, and clinical covariates. Among 21,447 patients included in the study, eosinophil count increased by a mean of 2.6 cells/L (SD 157.2 cells/L) and HbA1c decreased by a mean of 0.82% (SD 2.09%) during follow-up. Increases in eosinophil count were associated with reductions in HbA1c (parameter estimate -0.366 per 1,000 cells/L; 95% CI -0.526 to -0.206; P < 0.001), whereas baseline eosinophil count was not significantly associated with HbA1c changes. The association remained significant in a sensitivity analysis excluding systemic glucocorticoid users. Increases in eosinophil count were associated with improved glycemic control in adults with diabetes, whereas baseline eosinophil count was not. These findings extend prior cross-sectional observations and suggest that changes in eosinophils may represent a biomarker of glycemic improvement. Further prospective studies are needed to determine whether this relationship reflects underlying immune-metabolic mechanisms or treatment-related changes.

Rationale: Conventional CT-derived pulmonary vascular metrics quantify vascular pruning but do not assess the spatial heterogeneity of vascular architecture. In COPD, peripheral microvascular pruning coexists with compensatory central vessel dilation, paradoxically preserving total vessel volume while reducing effective vascular penetration into lung parenchyma. Objectives: To develop and evaluate Vascular Proximity Coverage (VPC), a distance-based CT biomarker quantifying the spatial efficiency of pulmonary vascular distribution. Methods: This retrospective analysis included 7,881 COPDGene participants [PHASE] with complete vascular segmentations. VPC was estimated by integrating the area under the cumulative lung tissue coverage curve at distances of 0-10 mm from segmented vessel voxels. Multivariable models assessed associations with lung function, six-minute walk distance (6MWD), St. George's Respiratory Questionnaire (SGRQ) score, mMRC dyspnea score, and all-cause mortality over a median follow-up of [10.8] years, adjusting for demographics, smoking, emphysema, scanner type, and BV5. The cohort was stratified at median BV5 and VPC into concordant (high/high, low/low) and discordant (high-BV5/low-VPC, low-BV5/high-VPC) vascular phenotypes. Results: VPC declined across disease severity (5.82 {+/-} 0.50 [UNITS] in non-smokers to 5.21 {+/-} 0.64 in GOLD 4; p < 0.001). Among spirometrically normal smokers, current smokers had lower VPC than former smokers (5.31 vs. 5.68; p < 0.001). Each standard deviation decrease in VPC was independently associated with reduced 6MWD ({beta} = -16.8 m, 95% CI: -19.5 to -14.0), worse SGRQ ({beta} = 1.9, 95% CI: 1.4 to 2.4), and greater dyspnea (RR = 1.05, 95% CI: 1.03 to 1.08; all p < 0.001). The lowest VPC quartile was associated with 38% higher mortality than the reference quartile (HR = 1.38, 95% CI: 1.18 to 1.62; p < 0.001). The high-BV5/low-VPC phenotype showed 24% higher mortality than high-BV5/high-VPC (HR = 1.24, 95% CI: 1.08 to 1.42; p = 0.002) despite preserved vessel volume. Conclusions: VPC captures spatial heterogeneity of vascular architecture independent of vessel volume, is reduced in smokers without airflow obstruction, and identifies excess mortality risk among individuals whose total vessel volume appears preserved.

Flow cytometry is central to clinical hematology, generating extensive archives comprising billions of single-cell profiles linked to expert-derived diagnostic summaries. While these records represent a largely untapped repository of collective expert knowledge, the underlying single-cell annotations required for supervised machine learning are typically not recorded in a standardized, accessible format. We present IMPGAT (implicit gating), a weakly supervised multiple-instance learning framework that reconstructs single-cell gating strategies directly from existing sample-level cell-type frequency distributions. The model utilizes an induced set attention mechanism to map single-cell marker profiles to specific cell populations, requiring neither explicit cell-level labels nor unsupervised clustering during training. We evaluate the framework on acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML) cohorts, totaling 711 measurements from 268 individuals and 240 million single cells. In ALL, the model's quantification of cell-type abundances shows no statistically significant difference from human inter-expert variability. In AML, the framework achieves blast quantification on par with a recently established supervised state-of-the-art model that relies on explicit single-cell annotations for training. Error analysis shows that mismatches between model predictions and expert gating primarily involve expert-labeled non-cellular events, reflecting the non-incorporation of any explicit physical prior knowledge (such as the geometric relationships traditionally used to identify doublets and debris in scatter cytograms) in model training. Conversely, model certainty correlates with phenotypic borderline cases where manual bisection is inherently difficult. Our approach enables the systematic utilization of retrospective clinical records for the development of automated, standardized gating logic and is applicable to broader high-throughput contexts, such as large-scale immune monitoring or phenotypic drug screening.

Background: Conventional discrimination metrics, such as the area under the receiver operating characteristic curve (AUROC), do not characterize the composition of outcomes among patients flagged as high risk by clinical prediction models. This limitation is particularly relevant in critical care, where competing deaths and incomplete follow-up can complicate the interpretation of prediction alerts. We developed the Balanced Risk Integrity Index (BRII), a stratum-level metric defined as the proportion of target events among target events, competing deaths, and early-censored observations in a flagged high-risk stratum. BRII is algebraically related to positive predictive value (PPV) and is reported alongside its component decomposition to characterize the outcome composition of flagged patients. Methods: We evaluated BRII using seven acute kidney injury (AKI) prediction models developed in MIMIC-IV and externally validated in MIMIC-III and eICU-CRD (116,010 patients in total). Additional simulations examined the effects of competing deaths and early censoring on BRII and conventional model performance metrics. Results: For 48-hour AKI prediction in the top 10% risk stratum, BRII ranged from 0.90 to 0.96 in MIMIC-III, compared with less than 0.26 across all models in eICU-CRD, where 61.5% of evaluated patients were censored before the prediction horizon. In the eICU-CRD top 10% strata for the same target, 37.0% to 50.2% of patients were censored early, 12.2% to 15.5% had AKI, and 3.4% to 8.1% died before AKI. BRII and PPV ordered at least one model pair differently in 34 of 48 cohort-outcome-top-k groups. Under simulated random censoring for 48-hour AKI prediction, BRII declined by up to 0.80 while binary-label AUROC changed by less than 0.02. Conclusions: BRII provided complementary information on high-risk stratum composition, yielding different model rankings from PPV and revealing substantial differences across external validation cohorts. Interpreted alongside PPV and its component decomposition, BRII offers a standardized approach to evaluating clinical prediction alerts across models and observation environments.

Asthma patients suffer from neuropsychiatric disorders (NPDs) more often than expected and NPD patients have a higher prevalence of smoking, which can lead to asthma. To investigate the genetic mechanisms underlying this complex relationship, we performed an integrative analysis of GWAS summary statistics from four asthma subtypes, 12 NPDs, and six substance use behaviors (SUBs). Asthma shares an important genetic component with post-traumatic stress disorder (PTSD), major depressive disorder, and attention-deficit/hyperactivity disorder (ADHD) that cannot be fully accounted for by smoking. Multi-trait GWAS identified 354 independent variants mapping to 178 pleiotropic genes capturing the core of the shared genetics of asthma, NPD, and SUB. Single-cell transcriptomic analysis across blood, brain, and lung tissues revealed pleiotropic gene overexpression in 57 cell types and three candidate mechanisms linking asthma and NPDs. First, immune-related pleiotropic genes are expressed in microglia, blood and lung immune cells, suggesting shared immune dysfunction. Second, the pleiotropic genes EP300, RERE, and MAML3 regulate NOTCH signaling and are expressed in neurons, microglia, and lymphoid progenitors. Third, pleiotropic genes expressed in ciliated cells of both lung and brain implicate ciliary dysfunction. These findings may help identify new treatment strategies for patients with co-occurring asthma and neuropsychiatric disorders.